Module 28: The Kepler Problem: Planetary Mechanics 

28.1 Introduction Kepler's Laws:^ 

1 . Each planet moves in an ellipse with the sun at one focus. 

2. The radius vector from the sun to a planet sweeps out equal areas in equal time. 

3. The period of revolution T of a planet about the sun is related to the major axis 
A of the ellipse by 

T^ = kA^ (28.1.1) 

where k is the same for all planets. 

28.2 Planetary Orbits: The Kepler Problem 

Since Johannes Kepler first formulated the laws that describe planetary motion, scientists 
endeavored to solve for the equation of motion of the planets. In his honor, this problem 
has been named The Kepler Problem. 

When there are more than two bodies, the problem becomes impossible to solve 
exactly. The most important "three-body problem" at the time involved finding the 
motion of the moon, since the moon interacts gravitationally with both the sun and the 
earth. Newton realized that if the exact position of the moon were known, the longitude 
of any observer on the earth could be determined by measuring the moon's position with 
respect to the stars. 

In the eighteenth century, Leonhard Euler and other mathematicians spent many 
years trying to solve the three-body problem, and they raised a deeper question. Do the 
small contributions from the gravitational interactions of all the planets make the 
planetary system unstable over long periods of time? At the end of 1 8th century, Pierre 
Simon Laplace and others found a series solution to this stability question, but it was 
unknown whether or not the series solution converged after a long period of time. Henri 
Poincare proved that the series actually diverged. 

Poincare went on to invent new mathematical methods that produced the modem 
fields of differential geometry and topology in order to answer the stability question 

using geometric arguments, rather than analytic methods. Poincare and others did manage 
to show that the three-body problem was indeed stable, due to the existence of periodic 
Solution. Just as in the time of Newton and Leibniz and the invention of calculus. 



^ As stated in An Introduction to Mechanics, Daniel Kleppner and Robert Kolenkow, McGraw-Hill, 1973, 
p401. 
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unsolved problems in celestial mechanics became the experimental laboratory for the 
discovery of new mathematics. 



28.3 Reducing the Two-Body Problem into a One-Body Problem 

We shall begin our solution of the two-body problem by showing how the motion of two 
bodies interacting via a gravitational force (two-body problem) is mathematically 
equivalent to the motion of a single body with a reduced mass given by 

^=^^ (28.3.1) 
+ 

that is acted on by an external central gravitational force. Once we solve for the motion of 
the reduced body in this equivalent one-body problem, we can then return to the real two- 
body problem and solve for the actual motion of the two original bodies. 

The reduced mass was introduced in Section 10.7 of these notes. That section used 
similar but different notation from that used in this chapter. 

Consider the gravitational force between two bodies with masses and as 
shown in Figure 28.1. 

body 1 
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body 2 

Figure 28.1 Gravitational force between two bodies. 

Choose a coordinate system with a choice of origin such that body 1 has position Fj and 
body 2 has position (Figure 28.2). The relative position vector r pointing from body 2 
to body 1 is r = Fj - . We denote the magnitude of r by |r| = r , where r is the distance 
between the bodies, and f is the unit vector pointing from body 2 to body 1, so that 

r = r f (28.3.2) 



2 




body 2 



Figure 28.2 Coordinate system for the two-body problem. 

The force on body 1 (due to the interaction of the two bodies) can be described as 



i,2=-^i,2r = -G^T^r, 



(28.3.3) 



Recall that Newton's Third Law requires that the force on body 2 is equal in magnitude 
and opposite in direction to the force on body 1, 



F = -F 

*1,2 *2,1' 



Newton's Second Law can be applied individually to the two bodies: 

2-: 



(28.3.4) 



F, , = m, — ^ 



(28.3.5) 
(28.3.6) 



Dividing through by the mass in each of Equations (28.3.5) and (28.3.6) yields 

Fu2 d'r, 



(28.3.7) 



^2 dt 
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(28.3.8) 



Subtracting the expression in Equation (28.3.8) from that in Equation (28.3.7) gives 
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d\ d% J^r 
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de df 



de 



(28.3.9) 



Using Newton's Third Law as given in Equation (28.3.4), Equation (28.3.9) becomes 



1,2 



— + — 

m, m 



2/ 



d^V 

de 



(28.3.10) 



Using the reduced mass jU , as defined in Equation (28.3.1) , 

- = — + —, (28.3.11) 

Equation (28.3.10) becomes 

F 

1,2 

F 

1,2 

where Fj 2 is given by Equation (28.3.3). 

Our result has a special interpretation using Newton's Second Law. Let be the 

reduced mass of a reduced body with position vector r = r f with respect to an origin O , 

where ris the unit vector pointing from the origin O to the reduced body. Then the 
equation of motion, Equation (28.3.12), implies that the body of reduced mass jU is under 

the influence of an atfractive gravitational force pointing toward the origin. So, the 
original two-body gravitational problem has now been reduced to an equivalent one-body 
problem, involving a reduced body with reduced mass 11 under the influence of a central 

force - Fj 2 f . Note that in this reformulation, there is no body located at the central point 

(the origin O ). The parameter r in the two-body problem is the relative distance 
between the original two bodies, while the same parameter r in the one-body problem is 
the distance between the reduced body and the central point. 



d^ 

de 

d^f 



dt 



(28.3.12) 



28.4 Energy and Angular Momentum, Constants of the Motion 

Consider the reduced body with reduced mass given by Equation (28.3.1), orbiting about 

a central point under the influence of a radially attractive force given by Equation 
(28.3.3). The equivalent one-body problem has two constants of the motion, energy E 
and the angular momentum L about the origin O . Energy is a constant because there are 
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no external forces acting on the reduced body, and angular momentum is constant about 
the origin because the only force is directed towards the origin, and hence the torque 
about the origin due to that force is zero (the vector from the origin to the reduced body 
is anti-parallel to the force vector and sinjr = ). Since angular momentum is constant, 
the orbit of the reduced body lies in a plane with the angular momentum vector pointing 
perpendicular to this plane. 

Choose polar coordinates (r, 6) for the reduced body (see Figure 28.3), where r is the 
distance of the reduced body from the central point that is now taken as the origin, 6 is 
the angle that the reduced body makes with respect to a chosen direction, and which 
increases positively in the counterclockwise direction, and F = - F^^r . 




Figure 28,3 Coordinate system for the orbit of the reduced body. 

Since the force is conservative, the potential energy with choice of zero reference point 
U (oo) = is given by 

Uir) = -^^. (28.4.1) 



The total energy E is constant, and the sum of the kinetic energy and the potential 
energy is 

£ = l^v^-^». (28.4.2) 
2 r 

The kinetic energy term, ^v^ 1 2 , has the reduced mass and the relative speed v of the 
two bodies. The velocity is cylindrical coordinates is given by (add link) 



V = V . f + V, 6, 

rad tan ' 



2 2 . 2 . ... 

V =v. =1 — I +ir 

lad tan ■ - ■ 



(dr_Yj^ddY (28.4.3) 
[dtj [ dtj 



where v^^^ = dr I dt and v^^ = r{dd I dt) . Equation (28.A.2) then becomes 
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E = —u 
2 



(dr_\ 
\dt} 



I 



dd\' 
dt } 



(28.4.4) 



The magnitude of the angular momentum with respect to the center of mass is 



dt 



(28.4.5) 



28.5 The Orbit Equation for the Reduced Body 

There are two approaches to describing the motion of the reduced body. We can 
try to find both the distance from the origin, r{t) and the angle, d{t) , as functions of the 
parameter time, but in most cases explicit functions can't be found analytically. We can 
also find the distance from the origin, r{d ) , as a function of the angle 6 . This second 
approach offers a spatial description of the motion of the reduced body. 



We begin with Newton's Second Law 



F = /^a 

( 

r = ij. 



d^r 



df 



- r 



ldd\'^. 

\dt} y 



(28.5.1) 



Note that the motion is no longer circular, the radius is not constant so there is a term in 
the radial acceleration d'^'r I dt^ along with the usual centripetal acceleration term 
-r{d6 1 dt f . Setting the components equal, using the constant of motion 
L = nr^{dd I dt) and rearranging, 



d^r 



L Gm^m^ 



dt fir 



(28.5.2) 



What we will do is use the same substitution u=l/r and change the independent 
variable from t to r , using the chain rule twice, since Equation (28.5.2) is a second-order 
equation. That is, the first time derivative is 



dr dr du dr du dd 
dt du dt du dd dt 



(28.5.3) 



From r = l/u we have drldu = -\lu^ . Combining with the angular velocity 

dd I dt = LI jur^ in terms of L and u , dd I dt = Lu^ I n , Equation (28.5.3) becomes 
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dr 1 du Lu^ 
dt dd n 

Taking the second derivative with respect to t , 



du L 
dd fi 
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d\L(L .\ 



de^ II 
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I 



d u 2 L 
dd^ 11^ 

Substituting into Equation (28.5.2), and using r = \lu , 



= — u - (jm^m^ u . 



-u 



dd n n 



Canceling the common factor of and rearranging, 



d^u uGm.m^ 
+ u = - 



dd' 



(28.5.4) 



(28.5.5) 



(28.5.6) 



(28.5.7) 



Equation (28.5.7) is mathematically equivalent to the Harmonic Oscillator Equation with 
a constant term. The solution consists of two parts: the angle-independent solution 



Wo = 



(28.5.8) 



and a sinusoidally varying term of the form 



(28.5.9) 



where A and 6^ are constants determined by the form of the orbit. The expression in 

Equation (28.5.8) is the inhomogeneous solution and represents a circular orbit. The 
expression in Equation (28.5.9) is the homogeneous solution (as hinted by the subscript) 
and must have two independent constants. Define 



(28.5.10) 



This is called the semilatus rectum. Then = 1 / and we can write the solution to 

Equation (28.5.7) as the sum of the inhomogeneous (Eq. (28.5.8)) and homogenous (Eq. 
(28.5.9)) pieces 
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" = "o + "h = -(1 + r,A(cos(d - d^))) 



r 

r= 



(28.5.11) 



We have two constants to choose. Define A = e I where 

e = ^1 + lixEr^ / =^l + 2EL^ I \x{Gm^ mj" (28.5.12) 

a quantity called eccentricity. Choose 0^=^ then cos(0 -7t) = - cos(0) . Thus Equation 
(28.5.1 1) can be written in the form 

r . (28.5.13) 

1 - ecos0 

The two constants of the motion in terms of and e are 

1 



E = 



2r 



The orbit equation as given in Equation (28.5.13) is a general conic section and is 
perhaps somewhat more familiar in Cartesian coordinates. Let x = rcosd and 
j; = rsin0 , with = + . The orbit equation can be rewritten as 

r = r^+e rcosd . (28.5.15) 
Using the Cartesian substitutions for x and y , rewrite Equation (28.5.15) as 

(x^ + yy =r^ + ex. (28.5.16) 
Squaring both sides of Equation (28.5.16), 

jc^+j;2 =rQ2 + 2exro+eV. (28.5.17) 

After rearranging terms, Equation (28.5.17) is the general expression of a conic section 
with axis on the x -axis, 

x'(l-e')-2exro + / = r^' (28.5.18) 
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(we now see that the dotted axis in Figure 28.3 can be taken to be the x -axis). 

For a given > , corresponding to a given nonzero angular momentum as in Equation 
(28.5.10), there are four cases determined by the value of the eccentricity. 

Case 1: When £ = , E = < and r = r^. Equation (28.5.18) is the equation for a 
circle, 

x'+/=ro' (28.5.19) 
Case 2: When < £ < 1 , E^^ < E <0 and Equation (28.5.18) describes an ellipse, 

y^ + Ax^-Bx = k (28.5.20) 

where A>0 and ^ is a positive constant. (Appendix 17.C shows how this expression 
may be expressed in the more traditional form involving the coordinates of the center of 
the ellipse and the semimajor and semiminor axes.) 

Case 3: When £ = 1 , E = and Equation (28.5.18) describes a parabola, 

x = ^ — ^. (28.5.21) 

Case 4: When £ > 1 , £" > and Equation (28.5.18) describes a hyperbola, 

y^-Ax^-Bx = k (28.5.22) 
where ^ > and A: is a positive constant. 

28.6 Energy Diagram, Effective Potential Energy, and Orbits 

The energy (Equation (28. A. 7)) of the reduced body moving in two dimensions can be 
reinterpreted as the energy of a reduced body moving in one dimension, the radial 
direction r , in an effective potential energy given by two terms, 

2jUr r 

The total energy is still the same, but our interpretation has changed; 



E = Kn+U^,,=^^ ^1 (28.6.2) 



^ dt ) Ipir 
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where the effective kinetic energy K^^^ associated with the one-dimensional motion is 



1 (dr 



2 \dt 



(28.6.3) 



The graph of U^^^ as a flinction of T = rlr^^, where as given in Equation (28.5.10), is 
shown in Figure 28.4. The upper curve (red, if you can see this in color) is proportional 
to l} I {2fxr^) ~ll 2r^ . The lower blue curve is proportional to -Gm^m^ I r — 1/r . The 
sum U^^f is represented by the green curve. The minimum value of 
will be shown analytically below. The vertical scale is in units of -U^^^ (r^ ) 



IS at r = , as 




Figure 28.4 Graph of effective potential energy. 

Whenever the one-dimensional kinetic energy is zero, K^^^ = , the energy is equal to the 
effective potential energy. 



(28.6.4) 



Recall that the potential energy is defined to be the negative integral of the work done by 
the force. For our reduction to a one -body problem, using the effective potential, we will 
introduce an effective force such that 



B 



B 



U,,,^,-U^^^,=-f¥^''-df = -fFf'dr 

A A 



(28.6.5) 
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The fundamental theorem of calculus (for one variable) then states that the integral of the 
derivative of the effective potential energy function between two points is the effective 
potential energy difference between those two points, 



\dU. 



U.rf,B-U,,,,A=J-^dr (28.6.6) 



Comparing Equation (28.6.6) to Equation (28.6.5) shows that the radial component of the 
effective force is the negative of the derivative of the effective potential energy, 

pcs^_dU^ (28.6.7) 
dr 

The effective potential energy describes the potential energy for a reduced body moving 
in one dimension. (Note that the effective potential energy is only a function of the 
variable r and is independent of the variable 6 ). There are two contributions to the 
effective potential energy, and the total radial component of the force is 

F-.4t/...4[-^L-^] (28.6.8) 
dr dr\2jj,r r j 

Thus there are two "forces" acting on the reduced body, 

" ^r, centifugal ^r, gravity ' (28.6.9) 

with an effective centrifugal force given by 

d ( \ 



F =- 

r, centrifugal 



^Ijir^ j lir^ 



(28.6.10) 



and the conventional gravitational force 



Gm, 



^.gravity =-^^^- (28.6.11) 

r 



With this nomenclature, let's review the four cases presented in Section 17.3. 



Case 1: Circular Orbit E = 
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The lowest energy state, E^^ , corresponds to the minimum of the effective potential 
energy, E^^ = {U^^)^^ . When this condition is satisfied the effective kinetic energy is 



zero since E = K^^^ + U^^^ . The condition 



(28.6.12) 



implies that the radial velocity is zero, so the distance r from the central point is a 
constant. This is the condition for a circular orbit. The condition for the minimum of the 
effective potential energy is 



dU^n 
dr 



G n\n\ 



(28.6.13) 



We can solve Equation (28.6.13) for r , 



Gm^m^ 



(28.6.14) 



reproducing Equation (28.5.10). 
Case 2: Elliptic Orbit E^^<E<0 

When K^^ = , the mechanical energy is equal to the effective potential energy, E = U^^ , 

as in Equation (28.6.4). Having dr/dt = corresponds to a point of closest or furthest 
approach as seen in Figure 28.4. This condition corresponds to the minimum and 
maximum values of r for an elliptic orbit. 



E = 



Gm^ 



2 fir' 



(28.6.15) 



Equation (28.6.15) is a quadratic equation for the distance r , 



2 Gm, «7, L 
r + — -r = 



2^E 



(28.6.16) 



with two roots 



Gm, m 



1 '"2 



' ( G m^m2^ ^ 



111 



IE 



IE I IpLE 
Equation (28.6.17) may be simplified somewhat as 



(28.6.17) 
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Gfflj 1712 



( i 



IE 



1 + 



1I}E 



1/2 X 



(28.6.18) 



Recall from Equation (28.5.12), the square root is the eccentricity e , given by 

e = ^1 + IfiEr^^ I =^\ + lElJ I \x{G m^f (28.6.19) 



Thus Equation (28.6.18) becomes 



Gm, m 



r = — 



l'"2 



IE 



(28.6.20) 



A little algebra shows that 



l} I lxGm^m2 



I 



1- 



1 + 

l} I iiGm^m 



2L^E \ 



liiGm^m^) 



-2L El jU(GOT[OT2) 



(28.6.21) 



IE 



Substituting the last expression in (28.6.21) into Equation (28.6.20) gives an expression 
for the points of closest and furthest approach, 



r = — ^(l±e) 



The minus sign corresponds to the distance of closest approach. 



r = r ■ = 

imn 



(28.6.22) 



l + £ 



(28.6.23) 



and the plus sign corresponds to the distance of furthest approach. 



r = r = 

max 



l-£ 



(28.6.24) 



Case 3: Parabolic Orbit £" = 
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The effective potential energy, as given in Equation (28.6.1), approaches zero (1/^^^ ~* 0) 
when the distance r approaches infinity (r -» oo ). Since the total energy is zero, when 
r ^ 00 the kinetic energy also approaches zero, K^^^ . This corresponds to a 
parabolic orbit (see Equation (28.5.21)). Recall that in order for a body to escape from a 
planet, the body must have a total energy E = (we set U^^ = at infinity). This escape 
velocity condition corresponds to a parabolic orbit. 

For a parabolic orbit, the body also has a distance of closest approach. This 
distance r^^ can be found from the condition 

E = U^^= ^ = . (28.6.25) 

2jUr r 

Solving Equation (28.6.25) for r yields 



r,.=ir^ \r,; (28.6.26) 

the fact that the minimum distance to the origin (the focus of a parabola) is half the 
semilatus rectum is a well-known property of a parabola. 

Case 4: Hyperbolic Orbit £" > 

When E >Q , in the limit as r ^ oo the kinetic energy is positive, K^^^ > . This 

corresponds to a hyperbolic orbit (see Equation (28.5.22)). The condition for closest 
approach is similar to Equation (28.6.15) except that the energy is now positive. This 
implies that there is only one positive solution to the quadratic Equation (28.6.16), the 
distance of closest approach for the hyperbolic orbit 

^h^=7^. (28.6.27) 
1 + e 

The constant is independent of the energy and fi-om Equation (28.5.12) as the energy 

of the reduced body increases, the eccentricity increases, and hence from Equation 
(28.6.27), the distance of closest approach gets smaller. 



28.7 Orbits of the Two Bodies 

The orbit of the reduced body can be circular, elliptical, parabohc or hyperbolic, 
depending on the values of the two constants of the motion, the angular momentum and 
the energy. Once we have the explicit solution (in this discussion, r{d)) for the reduced 
body, we can find the actual orbits of the two bodies. 
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Choose a coordinate system as we did for the reduction of the two-body problem 
(Figure 28.5). 




Figure 28.5 Center of mass coordinate system. 
The center of mass of the system is given by 



(28.7.1) 



Let ?/ be the vector from the center of mass to body 1 and r^' be the vector from the 
center of mass to body 2. Then, by the geometry in Figure 28.5, 



(28.7.2) 



and hence 



»'l=»-l-Rem=»'l — — 



nt- 



(28.7.3) 



A similar calculation shows that 




(28.7.4) 



Thus each body undergoes a motion about the center of mass in the same manner that the 
reduced body moves about the central point given by Equation (28.5.13). The only 
difference is that the distance from either body to the center of mass is shortened by a 
factor jU / OT; . When the orbit of the reduced body is an ellipse, then the orbits of the two 
bodies are also ellipses, as shown in Figure 28.6. 
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Figure 28.6 The elliptical motion of bodies under mutual gravitation. 

When one mass is much smaller than the other, for example m^«.mj, then the reduced 
mass is approximately the smaller mass, 

pL = — L^ = ^^ = Wi (28.7.5) 



The center of mass is located approximately at the position of the larger mass, body 2 of 
mass m2 . Thus body 1 moves according to 




(28.7.6) 



and body 2 is approximately stationary, 

r^ = — ^r-^r^6 (28.7.7) 

^2 
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28.8 Kepler's Laws 
Elliptic Orbit Law 

Each planet moves in an ellipse with the sun at one focus. 



When the energy is negative, E <Q , and according to Equation (28.5.12), 



e = 



1 + - 



2EL' V 



(28.8.1) 



and the eccentricity must fall within the range ^ £ < 1 . These orbits are either circles or 
ellipses. Note the elliptic orbit law is only valid if we assume that there is only one 

central force acting. We are ignoring the gravitational interactions due to all the other 
bodies in the universe, a necessary approximation for our analytic solution. 

Equal Area Law 

The radius vector from the sun to a planet sweeps out equal areas in equal time. 
Using analytic geometry, the sum of the areas of the triangles in Figure 28.7 is given by 

\, . {r ^e) \, . (r ^e) 

AA = -{rAd)r + ^ ^ ArM = -{rAd)r + ^ ^ Ar (28.8.2) 



in the limit of small Ad 




Ar 

Figure 28.7 Kepler's equal area law. 
The average rate of the change of area, AA , in time. At , is given by 

\(rAd)r (rAd)Ar 
AA = —- — + — ' 



At 



At 



(28.8.3) 
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In the limit as At ->0 , Ad —>0 this becomes 



^ = ir^ (28.8.4) 

dt 2 dt 

Recall that according to Equation (28.A.6) (reproduced below as Equation (28.8.5)), the 
angular momentum is related to the angular velocity dd/dt by 

^-^ (28.8.5) 

dt pLr 

and Equation (28.8.4) is then 

^ = A. (28.8.6) 

dt 2jU 

Since L and n are constants, the rate of change of area with respect to time is a 

constant. This is often familiarly referred to by the expression: equal areas are swept out 
in equal times (see Kepler's Laws at the beginning of this chapter). 

Period Law 

The period of revolution T of a planet about the sun is related to the major axis A of 
the ellipse by 



where k is the same for all planets. 

When Kepler stated his period law for planetary orbits based on observation, he only 
noted the dependence on the larger mass of the sun. Since the mass of the sun is much 
greater than the mass of the planets, his observation is an excellent approximation. 

Equation (28.8.6) can be rewritten in the form 

2n — = L. (28.8.7) 
dt 

Equation (28.8.7) can be integrated as 

Jl^dA^j'^Ldt (28.8.8) 

orbit 

where T is the period of the orbit. For an ellipse. 
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Area = J dA = jtab (28.8.9) 



orbit 



where a is the semimajor axis and b is the semiminor axis. Thus we have 



7 = ^^. (28.8.10) 



Squaring Equation (28.8.10) then yields 



/l_2 2 27,2 

= V • (28.8.11) 



In Appendix 17.B, the angular momentum is given in terms of the semimajor axis and the 
eccentricity by Equation (28.B.10). Substitution for the angular momentum into Equation 
(28.8.11) yields 



— (28.8.12) 



jU GOTj a(l - £ ) 

In Appendix 17.B, the semi-minor axis is given by Equation (28.B.22) which upon 
substitution into Equation (28.8.12) yields 

^,^4^rVV ^28.8.13) 

Using Equation (28.3.1) for reduced mass, the square of the period of the orbit is 
proportional to the semi-major axis cubed, 

T' , " (28.8.14) 
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28.9 Worked Examples 



28.9.1 Example Satellite Orbits 

A satellite of mass is initially in a circular orbit of radius around the earth. The 
earth has mass » and radius . Let G denote the universal gravitational 
constant. Express all your answers in terms of R^, m^, m^, G and as needed. 




a) Find an expression for the speed Vq of the satellite when it is in the circular orbit. 

b) Find an expression for the mechanical energy E^^ of the satellite when it is in the 
circular orbit. Take U(r) -^0 as r ^ co . 



As a result of an orbital maneuver the satellite trajectory is changed to an elliptical orbit. 
This is accomplished by firing a rocket for a short time interval thus increasing the 
tangential speed of the satellite. The apogee (farthest distance from earth) of the elliptical 
orbit is three times the closest approach (perigee), 

= \ = \ ■ 

(In the figure, the small red circle represents the earth.) 

c) Use conservation of energy and angular momentum for the elliptic orbit to find an 
expression for the speed of the satellite, Vp , immediately after the rocket has 

finished firing. 



Solution: 
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a) This preliminary part should be found directly from Newton's Second Law and the 
Universal Law of Gravitation. The magnitude of the acceleration for the circular orbit is 
vIItq, and so 



m, — = G-^ ^ 



"s ^2 

b) The total mechanical energy is the sum of the kinetic energy and the gravitational 
potential energy, 

„mjn^ 1 2 
En = -G— ^ + -OT.v, 



2 



1 Q m,m^ 



c) Since r^= r^^, and = Sr^ , the condition that angular momentum is constant 

TpVp = r^v^ becomes r^v^ = (\)v^ , so = Vp / 3 . The condition that the mechanical 
energy is constant then becomes, 

1 2 ^fn^m^ 1 2 ^^s^e 

2 s p „ 2 r 



1 ^v,^ = 



2 H 3 



2 

v„ = — G 



_Q!!V!h 



3 To 



= ^(3/2)GmJr,. 



As a simple check, note that Vp > Vg . As a further check, some minor algebra shows that 
after the rocket bum, the final mechanical energy is = -Gm^m^ I (4ro) = -Gm/n^ I A , 
where A = rQ+ 3rQ is the major axis of the ellipse. 
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28.9.2 Example Halley's Comet 

The equation for any orbit in an inverse square gravitational field is given by 



where 



r = ■ 



^0 



1 - e cos 6 



^0 = 



IxGm^m^ 



(28.9.1) 



(28.9.2) 




In Equation (28.9.2), L is the angular momentum, ^ = m^m2l{m^ + ^2) is the reduced 

mass, and e is the eccentricity of the orbit. When < £ < 1 , the orbit is an ellipse (in the 
above figure, e = 3/5). The period, T , depends only on the length of the major axis, A , 
of the ellipse, which is given by 



A = 



n 



(28.9.3) 



Halley's Comet is in an elliptic orbit about the sun. The eccentricity of the orbit is 
£ =0.967 and the period is T = 16y . The mass of the sun is m^^ = 1.99xl0'''kg . The 
mass of Halley's Comet is negligible compared to the sun. 

a) Using this data, determine the distance of Halley's Comet at closest approach 
(perihelion) to the sun, and furthest distance (aphelion) from the sun. 

b) What is the speed v of Halley's Comet when it is closest to the sun? 



Solution: 
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Before diving into the numerical calculations, let's do some preliminary math. 

First, note that when m^^^^^ «: m^^ (note that m^^^^^ is not given in the problem), 
m^+m^-^m^ and 

^ = = — . (28.9.4) 

«sun+ '"comet 1 + '"comet / "^sun 



Next, from Equation (28.9.1), we have 



'•p=7^, ';=T^ (28.9.5) 



1+e 1-e 
and combining with Equation (28.9.3) gives 



'•p=|'-o(l-0' ';=f'-o(l + 0- (28.9.6) 



As a quick check, note that r^+r^ = A, the major axis 



Next, anticipating part (b), we expect to use angular momentum considerations. At 
perihelion, the comet must be moving perpendicular to the vector from the sun to the 
comet, and so the magnitude of the angular momentum in terms of Vp , and m^^^ is 

^ = '"comet Vp. (28.9.7) 

Combining with Equation (28.9.2) (with =in^^^, =''^comet' vice versa) and the 
simplification for jU as given in Equation (28.9.4) and solving for Vp gives 

^ _ L _ ylGr^m^^ml^^^ _ V^'-Q'^sym _ / g(l + £K^ (28 9 8) 

'"comet'p '"comet'p 'p \ 'p 

The mass of the comet does indeed drop out of this problem if m^g^^^ «: m^^ . 

It's time to do the numbers. The calculations presented here were done by computer, 
keeping almost arbifrary precision in the intermediate calculations, and rounded to three 
figures (even though the period is given to only two figures). For the Newtonian 
gravitational constant, G = 6.67 x 10~" N • m^ • kg"^ was used. 
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a) Since the major axis A is used for determining both and , find that quantity first. 
From Equation (28.9.3), with ^1 + ^2= m^^^ , 



1/3 



'2(76yx3.16xl0's-y-')(6.67xl0-"N-m'-kg-')(l.99xl0'°kgy 



1/3 



\ 

= 5.37xl0"m 



(28.9.9) 



from which Equation (28.9.6) gives 



=8.86x10'° m, =5.28x10'' m. 



(28.9.10) 



These results are roughly half and thirty times the earth-sun distance, respectively; is 
roughly the distance from the sun to the ex-planet Pluto. If fact, the period of Halley's 
Comet is roughly 1/ Vs the period of Pluto's orbit, consistent with Equation (28.9.3). 



A graph of the orbit is shown here: 



The tiny dot (red, if viewed in color) represents the sun and is not to scale; a circle 
representing the sun to scale on this scale is too small to be seen (^p > lOOi?^^^^). 

b) Equation (28.9.8) then gives, with the result in Equation (28.9.10), 



v„ = 



G(1 + £K 



= 5.43xl0'm-s"' 



" A r ■ . (28.9.11) 



p 



This is essentially (but of course smaller than) the escape velocity from the sun. In fact 
it's not too hard to show that Vp / v^'^^^p^ = -^(l + £ ) , which is 0.9835 for the eccentricity of 
this orbit. 



24 



28,9,3 Example Lunar Orbit Collision 



A lunar mapping satellite of mass is in a circular orbit around the moon, and the orbit 
has radius = 1 .5 R^^ where is the radius of the moon. A repair robot of mass 
m,. < is injected into that orbit, but due to a NASA sign error it orbits in the opposite 

direction. The two collide and stick together in a useless metal mass. The point of this 
problem is to find whether they create more junk orbiting the moon or crash into the lunar 
surface. The mass of the moon is denoted by m„ . The universal gravitational constant is 
denoted by G. 




a. What is the initial orbital velocity of the mapping satellite, Vg ? Express your 
answer in terms of R^ , , and G. 



h. What is the speed of the space junk (satellite and robot) immediately after the 
collision? Write it as / , where you must determine the number / . Express 
your answer in terms of and m,. . 



c. After the collision, the orbit of the space junk has changed. Use conservation of 
energy and angular momentum to find an equation for the closest approach r^^ of 

the space junk to the moon. 

d. Solve your equation in part c for the number / assuming that the closest 
approach = R„, the space junk hits the moon. 



Solution: 



a) The speed of the mapping satellite undergoing uniform circular motion can be found 
from the force equation. 



Gm m mv, 

2 



2 

(28.9.12) 
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So the speed is 



\ = ^GmJR^ (28.9.13) 
b) The momentum flow diagram for the collision is shown below. 



Because there are no external forces, the momentum is constant and so 

"^/o ~ '^r^o = ('^. + ^r^^f ■ (28.9.14) 

Thus the speed after the collision is 

m -m 

V =^ ^v„ = /v„, (28.9.15) 

s r 

where the ratio of speed after the collision to the speed before the colhsion is given by the 
number 

/ = ^ = ^ ^. (28.9.16) 

s r 

c) After the collision, the energy equation is given by 

1 - G(m+m^m„ 1 . G(m+m^m„ 

-(m, + m^v/ ^ = -(m, + mjv' ^.(28.9.17) 



Setting Vf = /v^ and simplifying yields 



1 . . Gm„ 1 - Gm„ 

-f\ - = -v„^ (28.9.18) 

2-^ " /?„ 2 ^ r ^ ^ 



The angular momentum equation is 
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(28.9.19) 



Again setting v. = /v^, Eq. (28.9.19) becomes 



(28.9.20) 



Eq. (28.9.20) implies that = RqJvq I which we can substitute into Eq. (28.9.18) 
yielding 



1 /■2 2 



Gm^ \lR,jv,\ 



r 

V P I 



Collecting terms yields 



1 /■2 2 



1 -^0 

2 

\ P I 



- Gm^ 



- -\ 



(28.9.21) 



(28.9.22) 



If we assume that at the closest approach =R^, the space junk hits the moon. Then 
using the values =R^, v„ = ^Gm^ I , and R^ = (3/2)R^ , Eq. (28.9.22) becomes 



Gm^ 



3R 



M 



I 9R^\ 
AR ^ 



= Gm„ 



( 2 



1 \ 



V^^M Rm I 



(28.9.23) 



We can solve Eq. (28.9.23) for / : 



f = 4AI5 



(28.9.24) 
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28,9.4 Example Lowest Energy Solution 



The effective potential energy for a linear restoring central force F = -krr is given by 

•-^ effective ^ 2 ' 

2ixr 2 

Find the radius and the energy for the motion with the lowest energy. What type of 
motion does this correspond to? If the energy is slightly greater than the lowest energy, 
what type of motion would that correspond to? 
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Solution: 

Taking the derivative of the effective potential with respect to the radius r and setting the 
derivative equal to zero at , 



+ A: r„ = 



' '0 

, T , 1/4 



The minimum energy is 



„ . . 1} Ikji f- k( L' 



T ,1/2 



The orbit at this radius is a circle. 
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For a slightly larger energy, the orbit will oscillate about the minimum radius. Although 
it takes a bit more work to show, and is not part of this problem, the orbit will be an 
ellipse with the origin at the center (not at a focus, as in a Keplerian orbit). 
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28.9.5 Example Effective Potential Energy 



A system of two particles with a reduced mass interacts via an attractive central force 

F = -ar^ f , where r is the relative distance between the two particles. The magnitude of 
the angular momentum for the equivalent one particle problem about the origin is L . 

a) Find an expression for the effective potential energy and make a graph of the 
effective potential energy as a flinction of r. Find an expression for r that 
minimizes the effective potential energy. 

b) Indicate by a horizontal line on the graph sketch the total energy that would 

correspond to a constant relative distance r between the two particles. Find an 
expression for this energy. What type of motion does this correspond to for the 
equivalent one -particle problem? 



a) The potential energy associated with the given force, denoted U with no subscript, is 
given by 



where is the radius at which we choose to set the potential energy equal to zero. While 
we could choose any (finite) value for r^, for the purposes of making the graph we'll 
choose r^^Q . 

The effective potential energy is then 



From Equation (28.9.26) we see that the effective potential has both an attractive part 
(~ r"*) and the "repulsive" part ( — 1/r^). A plot of U^^^ as a function of r, with, for 
graphing purposes, a = L = ^ = 1 is shown below. Plots of U and the repulsive term are 
shown as well; U^^ is the upper curve, blue if viewed in color. 



Solution: 




(28.9.25) 



ZM a 4 ZM 

= —r + . 

2|Ur^ 4 2|Ur^ 



(28.9.26) 
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Differentiating Equation (28.9.26) and setting the derivative equal to zero at r = r^^ 
yields 



3 1 



= 



lib 



pia I 



In the above plot, with a = L = fj, = \, r^i„=l, consistent with the graph. 



b) Using the value of r^^-^^ from (28.9.27) in Equation (28.9.26) yields 



(28.9.27) 



2/3 



r2 / ^2 



1/3 



T 1/3 r4/3 

3 a L 



4 fx 



2/3 



(28.9.28) 



This value is shown as the black horizontal line in the plot below; with a = L = ^ = \ , 
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28.9.6 Example Transfer Orbit 

A space vehicle is in a circular orbit about the earth. The mass of the vehicle is and 
the radius of the orbit is 2R . It is desired to transfer the vehicle to a circular orbit of 

e 

radius AR . The mass of the earth is M . 




a) What is the difference in energy between the outer and inner circular orbits? 

b) An efficient way to accomplish the transfer is to use an elliptical orbit from point 
A at the inner circular orbit at to point B at the outer circular orbit (known as a 
Hohmann transfer orbit). This is accomplished by firing a rocket for a short time 
interval during each change of orbit thus increasing the tangential speed of the 
satellite. What velocity changes are required at the points of intersection, A 
and 5? 

c) Assume that the rocket bums fliel at a steady rate and the exhaust speed relative to 
the rocket is u . Using our early results for the rocket equation, how much fuel is 
burned at each of the rocket firings at the points of intersection, A and 5? 
Assume that before the firing at A , the mass of the rocket is . 
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Solution: 



a) Because « m^^ , the reduced mass n = m^. The total mechanical energy is the sum 
of the kinetic and potential energies, 

E = K + U 

1 2 ^m^M^ (28.9.29) 
= —my-G— — -. 

2 ' r 

For a circular orbit, the orbital speed and orbital radius must be related by the radial 
component of Newton's Second Law 

m M 

-G^-^ = -m — . (28.9.30) 



We can rewrite Eq. (28.9.30) as 



1 - 1 m M 

-mv^ = -G^ — (28.9.31) 

2 ^ 2 r ^ ^ 

Substituting Eq. (28.9.31) into Equation (28.9.29) gives 

E = _ = = lc/(r) . (28.9.32) 

2 r r 2 r 2 

In moving from a circular orbit of radius 2R^ to an orbit of radius AR^ , the total energy 
increases, as the energy becomes less negative. The change in energy is 



AE = E{r = 4R^)-E{r = 2Rj 



= _1g^-[-1g^1 (28.9.33) 
2 4i?. 2 2i? ' 



e 



8 R, 

b) The satellite must increase its speed at point A in order to move to the larger orbit 
radius and increase its speed again at point B to stay in a circular orbit. Denote the 
satellite speed at point A while in the circular orbit as . and after the speed increase (a 

"rocket bum") as ^ . Similarly, denote the satellite's speed when it first reaches point 

B as . and the speed of the satellite in the circular orbit at point B as ^ . The 

speeds . and v^^ are given by Equation (28.9.31). 



34 



While the satellite is moving from point A to point B (that is, during the transfer, after 
the first burn and before the second), both mechanical energy and angular momentum are 
conserved. Conservation of energy relates the speeds and radii by 



-f^S'^A f) -G— — - = -m(v„.) -G— — ^ 



(28.9.34) 



Conservation of angular momentum relates the speeds and radii by 

«s ^AJ (2^e) = ^B4 (4^e) ^ ^AJ = 2 . , 



(28.9.35) 



Substitution of Equation (28.9.35) into Equation (28.9.34) yields, after minor algebra. 



3 R 



'l GM^ 



(28.9.36) 



Equation (28.9.31) gives 



\\ GM^ 
2~Y 



'4 R 



(28.9.37) 



Thus the change in speeds at the respective points is given by 



/ \2 \\\ \GM^ 



and 



Av = V — V = \ a\— 



/ fl M GM 



(28.9.38) 



(28.9.39) 



Note that at both points, the speed must increase. 

c) Recall from our rocket equation (add correct link W08D2 Worked Example Rocket 
Problem ) that the change in the rocket speed is given by 



Av =v ^-v„ = Mln 

/- r,f rfi 



i 



(28.9.40) 



where the mass of the rocket after the burn is m = m„- Am^ with Am, > the amount 

of fiael that is necessary to bum in order to increase the speed by Av . We can write Eq. 
(28.9.40) as 
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Av.. . i^n- ^^f\ 

Exponentiating both sides yields 



= lni-^ ^i. (28.9.41) 



. , m.- Am, 
6^'''" = ^ /_ (28.9.42) 



We can now solve for hnij- 



Am =m„(l-e"^'''"). (28.9.43) 



7 

So for the first firing at point A we have that 

Am^^ = m„(l - e (28.9.44) 

where Av^ is given by Eq. (28.9.38). The amount of fuel after the firing is then 

m^^ = m^e (28.9.45) 
The amount of fuel burned during the second firing at point B is then 

^'^f,B = - ^~'^'''") = m^e-''^'\l - e'-''-) (28.9.46) 

where Av^ is given by Eq. (28.9.39). 
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Chapter 28 Appendices 

28.A: Derivation of the Orbit Equation 

Two presentations of the result of Equation (28.5.13) 

28.B: Properties of EUiptical Orbits 

The dynamical properties of objects in elliptical orbits 

28.C: Analytic Geometric Properties of Elhpses 

Demonstrating how the results of Appendices 17.A and 17.B are consistent with 
more familiar representations of ellipses 

28.D: Even More on Kepler Orbits 

Using geometry, vector algebra, but minimal calculus to find the orbit equation, 
introducing the Laplace-Runge-Lenz vector. 

Appendix 28.A: Derivation of the Orbit Equation Using 
Energy Methods 



Consider the reduced body with reduced mass given by Equation (28.3.1), orbiting about 
a central point under the influence of a radially attractive force given by Equation 
(28.3.3). Since the force is conservative, the potential energy with choice of zero 
reference point U (oo) = is given by 

U{r) = -^^. (28.A.1) 



The total energy E is constant, and the sum of the kinetic energy and the potential 
energy is 

E = -pLv^ (28.A.2) 

2 r 

The kinetic energy term, ^v^ 1 2 , has the reduced mass and the relative speed v of the 
two bodies. The velocity is cylindrical coordinates is given by (add link) 



V = V , f + V, 6, 

rad tan ' 



2 2 . 2 . ... 

V = v. +v =1 — I +ir 

rad tan ■ - ■ 



(dr_Y J^MY (28.A.3) 
[dtj [ dtj 



where v^.^^ = dr I dt and v^^^ = r{dd I dt) . Equation (28.A.2) then becomes 



37 



E = —u 
2 



(dr\ ( dd)' 
[dtj [ dtj 



(28.A.4) 



The magnitude of the angular momentum with respect to the center of mass is 



dt 



(28.A.5) 



We shall explicitly eliminate the 6 dependence from Equation (28.A.4) by using our 
expression in Equation (28.A.5), 



dd L 



dt iJ.r 



(28.A.6) 



The mechanical energy as expressed in Equation (28.A.4) then becomes 



Gm, OTo 



2 \dt ) 2 ixr 



(28.A.7) 



Equation (28. A. 7) is a separable differential equation involving the variable r as a 
function of time t and can be solved for the first derivative dridt. 



dt V jj, 



E -- 



1 ixr 



- + - 



(28.A.8) 



Instead of solving for the position of the reduced body as a function of time, we 
shall find a geometric description of the orbit by finding r{d) . We first divide Equation 

(28.A.6) by Equation (28.A.8) to obtain 



J0 

dr dr^ 
dt 



pir 



E- 



1 L Gm,m.Y 

— 2^ \ 

2 pir r I 



(28.A.9) 



The variables r and 6 are separable; 



dd = 



L 



i\lr') 



\ 2iur' ^ r ] 



1/2 



dr . 



(28. A. 10) 
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What follows involves a good deal of hindsight, allowing selection of convenient 
substitutions in the math in order to get a clean result. First, note the many factors of the 
reciprocal of r . Make the substitution u =l/r , du = r^) dr , with the result 



de = —^- ^. (28.A.11) 

t u +Gm^m2U 

[ I 

Experience in evaluating integrals suggests that we make the absolute value of the factor 
multiplying inside the square root equal to unity. That is, multiplying numerator and 
denominator by -sjl^/L, 

dd — — -. (28.A.12) 

As both a check and a motivation for the next steps, note that the left side {dd) of 
Equation (28.A.12) is dimensionless, and so the right side must be. This means that the 
factor of juGwijOTj/Z^ in the square root must have the same dimensions as m, or 

length"' ; so, define 

r^ = I fxGm^m^. (28.A.13) 

a quantity we previously encountered in Eq. (28.5.10) called the semilatus rectum. The 
differential equation then becomes 

de = - 



{2HE I U -u" +2ul rj 
du 



" 0^ (28.A.14) 

du 



{lliE I L + \l r^' -{u-\l r^Y) 



2\l/2 



du 



illiEr^'lL' + \-ir^u-\ff' 



Next, we note that the combination of terms l^Er^ 1 1} + \ is dimensionless, so define 

e = ^1 + luEr^"- 1 1} (28. A. 1 5) 
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a quantity we also previously encountered (Eq. (28.5.12)) called the eccentricity. The last 
expression in is then 



r, du 

de = - . / (28.A.16) 

(£ -{r^u-l)) 

From here, we'll combine a few calculus steps, going immediately to the substitution 
Tq M - 1 = £ cosa , du = -s sina da , yielding 

-e sincf da , ^ , _ 

^0 = - . , , , . (28.A.17) 

(e -e cos a) 

We can now integrate and find that 

d = Jdd = Jda = a + constant (28. A. 1 8) 



We have a choice in selecting the constant, and if we pick d=a-Jt , a=d+Jt , 
cosa = -cos 6 , the result is 

r = - = ^ (28.A.19) 

u \-scosd 

agreeing with our result in Equation (28.5.13). 

Note that if we chose the constant of integration to be zero, the result would be 

r = - = ^-^ (28.A.20) 

U 1 + £COS0 

which is the same trajectory reflected about the "vertical" axis in Figure 28.3, indeed the 
same as rotating by n . 



Appendix 28.B: Properties of an Elliptical Orbit 

We now consider the special case of an elliptical orbit. 
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Figure 28.B.1 Ellipse. 



In Figure 28.B.1, let a denote the semimajor axis, b denote the semiminor axis and 
denote the distance from the center of the ellipse to the origin of our coordinate system 
{r,d) . We shall now express the parameters a , b and in terms of the constants of the 

motion L , E , ^ , and . 
The semimajor axis: 

See Equation (28.A.19) above. The major axis A = 2a is given by 

^ = 2a = r_ + r„,„ (28.B.1) 
where the distance of furthest approach occurs when = 0, hence 

r^^^ = r(d = 0) = (28.B.2) 

and the distance of nearest approach occurs when 6 =jt , hence 

r . =r{d = it) = -^. (28.B.3) 

l+£ 



y 




















' ' / 

''max 




mm 





Figure 28.B.2: nearest and flirthest approach 



Thus 
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a = -(-^ + -^] = -^. (28.B.4) 

The semilatus rectum can be re-expressed in terms of the semimajor axis and the 
eccentricity, 

r„ = a(l-£'). (28.B.5) 

We can now express the distance of nearest approach, Equation (28. B. 3), in terms of the 
semimajor axis and the eccentricity. 



1 + £ 1 + £ 



r.=^ = — ^ = a(l-£). (28.B.6) 



In a similar fashion the distance of flirthest approach is 

^0 



1 - £ 1 - £ 



= a(l + £). (28.B.7) 



Figure 17.B.2 shows the distances of nearest and furthest approach. Using our resuhs for 
and £ from Equations (17.3.12) and (17.3.13), we have for the semimajor axis 



L' 1 

a = 



^Gm^m^^l_^l^2EL' I ix[Gm^m^ )')) (28.B.8) 
IE 

and so the energy is determined by the semimajor axis, 

£ = -^^. (28.B.9) 
2a 

The angular momentum is related to the semilatus rectum r^ by Equation (17.3.12). 
Using Equation (28. B. 5) for r^ , we can express the angular momentum (28. B. 3) in terms 
of the semimajor axis and the eccentricity, 

L = ^JfiGm^ m^r^ = yjfJ,Gm^ m^a{\ - e^) . (28.B.10) 

Note that 
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/ ^ (28.B.11) 



Location Xg of the center of the ellipse: 

From Figure 28.B.1 , the distance from a focal point to the center of the ellipse is 

Xo = r_-fl. (28.B.12) 

Using Equation (28. B. 7) for r^^^ , we have that 

Xg = + £)-£ = £ a. (28.B.13) 

Thus, from Equations (17.3.12), (28.B.13) and (28.B.8), 



Gm m 

x„ = £ a = - 







2E 



y]{\ + 2EL^ / ^{Gm^ m^f). (28.B.14) 



The semi-minor axis: 

From Figure 28.B.1, 



where 



which can be rewritten as 



Note that from Figure 28.B.1 , 



so that 



b=yl(r,'-x,') (28.B.15) 



r, = ^ , (28.B.16) 

1-£COS0^ 



-r^£cos0^ = Tg. (28.B.17) 



Xo=r,cos0,, (28.B.18) 



r^=ro + £Xo. (28.B.19) 
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Substituting Equation (28.B.13) for and Equation (28. B. 5) for into 

Equation (28.B.19) yields 



r = a(l - e^) + ae^ = a 



(28.B.20) 



The fact that = a is a well-known property of an ellipse reflected in the geometric 

construction, that the sum of the distances from the two foci to any point on the ellipse is 
a constant. Thus the semi-minor axis b becomes 

b = ^r/ - = Va' - e^a^ = ayjl - . (28.B.21) 



Using Equation (28. B. 11) for Vl - , we have for the semi-minor axis 



b = ^al} I pLGm^ 



(28.B.22) 



We can now use Equation (28. B. 8) for a in the above expression, yielding 



'^^'^^ I ^Gm^m^ =L (28.B.23) 



IE 



Speeds at nearest and furthest distances: 

At nearest approach the velocity vector is tangent to the orbit, so the angular 
momentum is 



L = ur-v„ 

^ mm p 



(28.B.24) 



and the speed at nearest approach is 



(28.B.25) 




Figure 28.B.3 Speeds at nearest and furthest approach 
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Using Equation (28.B.10) for the angular momentum and Equation (28.B.6) for 
Equation (28.B.25) becomes 



^ L _ ^iiGm^m^{\-e') _ \ Gm^m^(\-e') _ j^^^^^ B 26) 

A similar calculation show that the speed at furthest approach, 

^ L _ ^nGm^m^{\-e') _ i Gm^m^jl-e') _ i Gm^m^jl-e) ^^^^^^^ 

Appendix 28.C: Analytic Geometric Properties of Ellipses 

Consider Equation (28.5.18), and for now take e < 1 , so that the equation is that of an 
ellipse. It takes some, but not a great deal, of algebra to put this into the more 
familiar form 



b 



2 +Tr = l (28.C.1) 



where the ellipse has axes parallel to the x and y coordinate axes, center at (Xo,0), 
semimajor axis a and semiminor axis b . 

We rewrite Equation (28.5.18) as 

1-e 1-e 1-e 

We next complete the square. 



X + + ^ = — ^ + 



1-e' (1-e')' 1-e' 1-e' (1-e')' 

1 2 2 

x-^\ ^ y - 



(28.C.3) 



l-e'/l 1-e' (1-e')' 



We now multiply each side by the factor (1 - e')' / yielding 



= l (2B.C.4, 

1-8 / To (1 - e ) 



Eq. (28.C.4) simplifies to 
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(28.C.5) 



The last expression in (28.C.3) is the equation of an ellipse 



(x - eaf 



1 



(28.C.6) 



with semimajor axis 







(28.C.7) 



a = 



2 ' 



1-e 



semiminor axis 



b = 




(28.C.8) 



and center at (er^ / (1 - e^), 0) = (ea, 0) , as found in Equation (28.B. 13). 

Appendix 28.D: Even More on Kepler Orbits 

We've seen so far that for a Kepler Orbit, we have two constants of the motion: the 
angular momentum and the total energy. Since the angular momentum is a vector with 
three components, these constitute a total of four scalar constants of the motion. The 
Kepler Problem has six "degrees of freedom" (three position, three velocity), and so we 
expect to be able to find two more scalar constants of the motion. 

We might expect to be able to further identify any orbit by a vector in the plane of the 
orbit, perpendicular to the angular momentum, and this is indeed the case. Symmetry 
suggests that this vector would be along the major axis, and we'll see that this is the case 
as well. What follows uses a good deal of vector algebra, but minimal calculus, and leads 
to the orbit equation in a surprisingly simple form. 

We'll need two results from vector algebra that we haven't had to use yet. Specifically, 
for vectors a , b and c , we have 



These relations are not hard to derive in Cartesian coordinates; the derivations will not be 
reproduced here. As a check, however, note that the vector on the right side of the first 



a X (b X c) = b(a • c) - c(a • b) 
a • (b X c) = b • (c X a) = c • (a X b). 



(28.D.1) 
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relation is perpendicular to both a and b x c . If a , b and c are non-coplanar, the 
common magnitude of the scalars in the second relation is the expression for the volume 
of a parallelepiped with the three vectors forming the sides. 

Let's start with the known constant angular momentum, 

L=rxp = ;urxv (28.D.2) 

and re-express this quantity in a way that will allow us to use Newton's Laws. 
Specifically, consider the velocity in terms of polar coordinates, 

v = ^r + r— (28.D.3) 
dt dt 

so that the angular momentum can be expressed as 

L = ^(.r)x^.-^=^. ^rx-^ (28.D.4) 



fi'om which 

Jr _ J_L 

dt jU 



fx— = — -. (28.D.5) 



The advantage to this operation is that we now have an explicit scalar factor of \lr^ , 
which can and should be related to the same factor that appears in Newton's Law of 
Gravitation. However, in order to use that law, we need a vector relation involving f , 
and so we'll cross r into both sides of Equation (28.D.5), yielding 



rx 



dY\ r L 

— = — 
dt j r fi 



(28.D.6) 

dr \(d'\-- 
— = — r xL. 

dt GotjOTj \ dt j 



In the above, the first relation in Equation (28.D.1) was used to simplify the left side, and 
Newton's Law of Gravitation, in the form fx{d^r I dt^) = -{Gm^m^ I r^)r was used on the 
right side. Note the cancellation of the factor of the reduced mass n . 

Equation (28. D. 6) may now be integrated to obtain 
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r = 



- + A 



A = r- 



vxL 



(28.D.7) 



where A is a constant vector. Since r is in the plane of the orbit, and L is 
perpendicular to the plane of the orbit, A must lie in the plane of the orbit, as indicated 
above. Further, by considering extreme points of the orbit, where v 1 f , and hence 
vxL||r, we see that at these points A is in the direction parallel to the major axis. 
Since A is a constant vector, A must always be in this direction. By considering the 
vector A at perihelion (or at any point on the orbit of a circular orbit), we can see that the 

direction of A is that from the perihehon point to the focus; we'll need this result below, 
when we find the orbit equation. 

The magnitude of A is readily found by calculating 



A^ = A-A = r-r — f • (v x L) + ^- (v x L) • (v x L) .(28.D.8) 



Gm^m^ 



The first dot product is manifestly 1 . The middle term, the "cross term," is found using 
the second relation in Equation (28.D.1), 



- , 1 - 1 , 

r ■ (v X L) = L • (r X v) = — L • (r x = — L , 



(28.D.9) 



The third term is most easily evaluated by recalling that v ± L , so that vxL =vL and 
(v X L) • (v X L) = v^l} . Combining, we see that 



1 



1 



A' = \-— -L' + 



= 1 + 



2Z' (I 2 Gm^m^\ 



= 1 + 



12 

2L^E 



= e 



(28.D.10) 



Thus, the constant vector A is directed along the major axis and has magnitude equal to 
the eccentricity. 



The orbit equation is now found algebraically by taking the dot product of A and r : 
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A • r = COS0 = f • r r • (v x L) . (28.D. 1 1) 

The first term is merely the magnitude r of the position vector. The second term, 
repeating the calculation of Equation (28. D. 9) with r instead of f , is 

^- L^=r^, (28.D.12) 

with the result 

Arcosd =er cose =r-rQ. (28.D.13) 
Solving for r gives the orbit equation in the form 



°- . (28.D.14) 

1-£COS0 



It should be noted that what we call the vector A is a negative scalar multiple of the 
"Laplace-Runge-Lenz" vector (yes, it took three people to come up with this). 
Specifically, the L-R -L vector is in many sources given as 



-^LRL = P X L - Gm^m2 — = -/^ Gm^m^ A 

r 



(28.D.15) 



Our choice of the form for A allows 



= e and the direction of A to lead to 



Equation (28.D.11) without introduction of extra minus signs. 
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